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1 Introduction 


In this article we study 0{N) linear sigma models near criticality, within the 
framework of the Exact Renormalization Group (ERG, hereafter), mainly in 
its simpler nonperturbative approximation, the Local Potential Approximation 
(LPA). We concentrate mostly in the Wegner-Houghton equation although 
other equations are also considered [||, ||. 

Recent results @] have raised interesting questions about the behaviour of 
the RG near first-order phase transitions. They are in some conflict with exact 
results of Wegner-Houghton equation at ^ oo ||]. A systematic study of 
this subject at finite N is still lacking. Although our interest points in that 
direction, many aspects of second order phase transitions had to be previously 
worked out. This paper is meant to fill this gap, as well as review the derivation 
and projection of the equation with a modern language, paying special attention 
to some aspects not fully treated in the literature, which we hope will be use¬ 
ful for a subsequent study of 0{N) models at discontinuous phase transitions. 
Nevertheless, as this paper is focused in continuous transitions, we leave to the 
conclusions further comments on first-order ones. 

ERG methods are exact formulations of the RG in differential form (see 
Ref. 1^), that is, the evolution of the renormalized action along the RG flow 
is studied in terms of differential equations. There are many such equations, 
being one of the oldest that of Wegner and Houghton Q , which corresponds to 
a decimation in momentum space, i. e., a lowering of a sharp momentum cutoff. 
Unfortunately, those equations are very involved and exact solutions are only 
available for very simple models. Nevertheless, several feasible nonperturbative 
approximations are possible, being the most promising an expansion in powers 
of momenta. The LPA 0 is just the first term of such an expansion [^. 

The organization of the paper is as follows. In Section Wegner-Houghton 
equation is derived in some detail. The LPA is discussed in Section ^ and exten¬ 
sively worked out in Section where we compute critical indices of 0{N) models 
at finite N as well as its eigenoperators. It is also shown that irrelevant eigen¬ 
values, which are related to the anomalous dimensions of higher-dimensional 
composite operators, may be obtained with ease. The comparison with our re¬ 
sults from Polchinski equation as well as other determinations is left to Section 
H, while the special, exactly solvable, N ^ oo case is studied in Section^, where 
some striking results are obtained. Section ^ is devoted to the conclusions. Ei- 
nally, some peripheral subjects are treated in two Appendices. 

2 Derivation 

In this section a careful derivation of Wegner-Houghton equation is presented, 
trying to be as self-contained as possible. Nevertheless, for the sake of simplicity, 
we concentrate on the iV = 1 case. Once this case is mastered, its generalization 
for arbitrary N is straightforward. 

We assume that our action is regulated with a t-dependent sharp cutoff in 
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momentum space At = e“*Ao, with all modes with |g| > At already eliminated, 
where t parametrizes the RG flow and Aq is a fixed scale. A RG transformation 
consists of two steps, a blocking or elimination of short-distance degrees of 
freedom, and a change of length scale. In this way one goes from an action 
cutoff at At to one at At+r, in physical units. From now on, we work with 
dimensionless quantities, the dimensions being given, if needed, by At. We 
derive Wegner-Houghton equation first considering the blocking, and afterwards 
the rescaling. 

Let us parameterize our action as 


— 'y ^ / ''I ■ ■ ■ T<lm)4‘qi ' ' ' + 

m=2,4,... •^91 "'9"* 

V2{q)(l>q(l)-q + S[cj)], ( 1 ) 

where and 5q = (27r)'^(5(g). We have assumed Z 2 symmetry. The 

field is split as 

</., = 4°)+</>«, (2) 

where are the modes with |g| < e~'^ and the ones with e~'^ < jgl < 1. 
The blocked action is then defined through 





,-Sl4>] 


( 3 ) 


As we are seeking a differential equation for the renormalized action, we 
need to consider just an infinitesimal blocking, hence we discard terms of order 
or higher. For that purpose, it is convenient to rewrite the action, Eq. as 

^[<^] = sm + i f ■■■ f (4) 

^ ''q n=l "" "''JI 


with 


qM = _ 




...S(j3q 




( 5 ) 


and primes in integrals meaning momenta restricted to the shell e~'^ < |g| < 1. 

Written in the above form, the action is suited for treating the -integral 
in a diagrammatic expansion. Its Feynman rules are: 


1. Any n-legs vertex contributes with ''' fq ^qil--,qm where every leg is 
labelled by a momentum qi. 


2. Any propagator between a leg labelled by qi and another one labelled by 
92 is 5q^+q^V2^{qi). 

3. All gi-integrals are understood to affect the whole diagram and not only 
the contribution of its vertex. That is, every delta function always simpli¬ 
fies an integral. 
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4. Add the usual symmetry factors. 

Before going on we state two technical results. 

Lemma 1. For any function f{q;P), analytic around jgl, |P| = 1, it is 


J[dP]J^f{q;P)6,+p = 0{T), (6) 

where P stands for a finite number of momenta pi, with \pi\ < 1 , and J[dP] for 
integrals over all p^’s. 

Proof. Just integrate first over all P integrals. QED 

Lemma 2. For any function f{qi, q 2 ', P), analytic for any value of P and at 
kil, I 92 I = 1, it is 

[ Sq,+q^+pf{qi,q2;P)=0{T), (7) 

provided that P = 0 + 0{t). 

Proof. In spherical coordinates, integrating with respect to \q 2 \, 


/ / Sq^+q^+pf{qi,q2;P) 

Jqi Jq2 

^ f d\qi ' 

Je-^<\qi\<l (27r) 


911 


dLli 


( 8 ) 


dn 2 S{n)f{\qi\qi,\qi+P\q 2 -,P), 


with \qi + P\ being in the shell and 6(id) standing for the delta functions that 
constrain angular variables. The integral over |qi| will give a contribution of 
0{t). To convince ourselves that the rest of the integral is 0{t^) for P = 0, 
let us take 9i to be the angle between P and In this case, the angular 
contribution of Eq. ^ is 


/ 


dfli 


d^2 S{Ll) 



dOi 


d6»2 - 4 * 2 ) 


X J(|gi| sin0i - \qi + P| su\92)5{P + |gi| cos^i + \qi + P| cos02) (9) 


where <I>i is a short-hand for all the remaining angular variables and J is the 
appropriate Jacobian. From the last two delta functions, recalling that 1 — r < 

hi < 1 , 

01 =^- 02 +O(t), P = 0 + O(r). (10) 

QED 

Let us now compute the leading contribution, of a tree diagram 

assuming that each n-vertex is with [z] being the lowest integer greater 

or equal than z, propagators being 0{\/t) because of the delta function. Then, 

D — n{l2n+l2n-l) — P — n-{d2n+l2n-l)—nin = - l2n-l, (H) 

n n n n 

where is the number of n-vertices and P the number of propagators. Let us 
now put this result in a rigorous basis. 
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Figure 1: A general diagram consisting on n 1-vertices and 1 n-vertex. 


Lemma 3. The leading contribution in r of a tree-level diagram, 0{t^) is 

^ = 5 J2n^2n-l- 

Proof. Consider first a connected tree diagram consisting of one n-vertex 
and n 1-vertices, as in Fig. |l|. Its value is 



( 12 ) 


We perform the gi-integrals in pairs, using Lemma 2 to extract the leading 
dependences in r. At the end, we obtain, in agreement with Lemma 3, that it 
is for n even or else for n odd, with the final aid of Lemma 

1. (Recall the delta function hidden in Eq. H). 

Furthermore, let us assume that Lemma 3 is true for all diagrams with K 
vertices other than 1-leg ones. Let us now prove that it is also true for diagrams 
with K + 1 such vertices. To this end, note that a generic tree diagram Q with 
K + 1 more-than-one-leg vertices contains at least one n-vertex (for some n 7 ^ 1) 
with all but one leg connected to l-vertices.|^ 

If the selected vertex is a n-vertex, with n odd, perform all integrals asso¬ 
ciated with its legs but the one that connects it to the rest of the diagram. 
Performing the integrals in pairs, and using Lemma 2, it brings about a contri¬ 
bution It remains a diagram Q' with K more-than-one-leg vertices, 

where the former selected vertex enters now as an additional I-vertex (with 5*' 
replaced by some other function). See Fig. ||. 

^ To pick up one, just eliminate all 1-vertices, and take one 1-vertex of the remaining 
diagram (which must have at least one since it is again a tree). 
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nodd neven 

Figure 2: A vertex with n integrated legs, as explained in the proof of Lemma 
3. we have to distinguish between the n odd and the n even case. 


If n is even, then integrate over all its legs connected to 1-vertices but one. 
This is, using Lemma 2 again, There is left over a 2-vertex con¬ 

nected both to the remaining sub-diagram and to a 1-vertex. It can easily be 
seen that it is of the same order in r as one with the 1-vertex connected directly 
to the sub-diagram, without the 2-vertex. See again Fig. |^. 

Therefore, 


D{g) = D{g') + 


(13) 


and, by the induction hypothesis. 


D{g') 


1 

2 


E4-1 


i>l 

E + ^odd 


(14) 


Taking into account that /f = /f — (n — 1) -I- 1 and I^' = — 1, a proof of 

Lemma 3 follows. 

QED 

Note that every tree diagram must necessarily beQ Ii > 2, and, from Lemma 
3, it must be /i = 2 at leading order in r. Furthermore, in such case, any 
inclusion of a four-legs or higher vertex would imply at least one loop, and thus 
the only vertices allowed, apart from one-leg ones, are an arbitrary number of 
two-legs ones. Even more, these additional vertices, in order to be 0{t), are 
{Lemma 2) 



The argument that leads to the proof of Lemma 3 has to be slightly modified 
for diagrams with L loops. Any IPI part contains one g-integral per loop, 
with the only restriction that of being within the shell. After performing those 
integrals, we end up with a contribution 0{t^) times an effective tree diagram. 
Therefore, the correct formula is to add the number of loops to the result of 

^ Recall that every diagram must contain at least one propagator. 
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Lemma 3, but with Eq. computed considering any IPI part as a tree-level 
vertex. 

We conclude then that there cannot be diagrams with more than one loop 
and that one-loop diagrams must have zero external legs. Furthermore, it hap¬ 
pens again that diagrams with more than two legs generate further loops and 
that two-legs ones must have the form of Eq. 

All those considerations show that, at leading order in t, Eq. ^ is equivalent 
to ^ ^ 

5[</.] = + ^J' (16) 

The path integral, Eq. ||, is now Gaussian with the result 

+ K-g) + const, (17) 

where we assume that is positive definite and “const” means a possible 

term independent of . 

The other step in any RG transformation is a change of scale in all linear 
dimensions. This amounts to a standard dilatation transformation. 


k = e^q, 


. d+2-7, 


= e 




(18) 


which gives. 




Ul ^1 ) • • • 

X 

+ T fdS' - -— 


6 '^m(6 , 6 ^m) 

ci+2 —Tj 




(19) 

0 (r 2 ), 


where the prime in the last partial derivative merely indicates that it does not 
affect delta functions. 

Therefore, Wegner-Houghton equation is 


S=^ = \im 

dt T^o T 


= lim — 
T^o 2r 


+ dS- 


Uk 
2 - d - 77 


f s'_kSUsik,k) 

Jk 


-1 




■ , , d' SS 
VkK ■ wyT— -I- const. 
; uk d(pk 


( 20 ) 


To deal with the N arbitrary case, we define the partition function as 


Z[J*] = (21) 
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so that couplings are all finite at N —>■ oo. The derivation is completely analo¬ 
gous and the ERG equation reads 
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lim — 
r—>0 2r 


+ dS- 


1 

N __ 

2- d-r] 




SikSk 


i ( 

k \^-k,k) 


-1 


■ , (55 A ,, d' 6S 


( 22 ) 


where the trace is over flavour indices and 


C2l.---,2n 


5'^S 


. . . S(j)q 


0(1)=O 


(23) 


3 Projection 


The Local Potential Approximation (LPA) amounts to projecting the full action 
into a fixed kinetic term plus a general potential involving only the zero modes, 


5'=^y k'^'pk(l>-k + V{y)So, 


(24) 


where y = 4’h4'h- The LPA is, thus, just the first term in an expansion in powers 
of momenta |^. 

To project over constant fields, we define the operator V{x^) 0 which, acting 
on an arbitrary functional G, is 


rix^)G= exp (a:*— ) G 


(25) 


k=0 


The projection of the third term in the second line of Eq. ^ gives zero because 
of the explicit factor of k there. The first two terms in that line are simply 

dV{y)So + - —^— -x"^^^6o = [dV{y) + (2 - d - r])yV'{y)]So, (26) 

with x'^x^ = y. The second term of the right hand side of Eq. can be seen to 
cancel whereas the first one is 

p/ 


lim 


1 


-0 2Nt 


tr In 




The LPA amounts to 




S^^k^ + 


d'^V 

dx^dx^ 


So- 


(27) 


(28) 


Eq. 


I becomes, thus 
■A-d 


— ^ 


tr ln(l -I- m) -I- tr In ( -|- 2x''x^ j 


-I- u 


-b const 


= — [(A^ — 1) ln(l -b m) -b ln(l -b u -b ‘2yu')] + const. 


(29) 
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with u{y) = 2V'{y), 


_ 27T‘^/yTidl2) 

~ {27Ty 


(30) 


and r(z) the Euler Gamma function. 

The equation is further simplified if we perform one more derivative with 
respect to y, 


Soil 


Ad 


3u' + 2yu'' 
1 + M + 2yu 


- + {N-1) 


1 + u 


+ <5o[(2 — ri)u + (2 


d-ri)yu']. (31) 


If, instead, the projection is over the kinetic term and its normalization kept 
fixed under the RG flow, 

0 = -/ + + i(2 -d-r])p^, (32) 


it follows ry = 0. 

Moreover, the delta function in Eq. ^just reflects the difficulties of selecting 
one mode out of a continuum set. This problem is not present at finite volume, 
where the zero mode is neatly separated. In such a case the delta function 6o 
is just a coefficient proportional to the inverse number of modes. One can then 
rescale y y/So and Eq. ^ becomes 


Ad 

N 


+ (JV - 1 )^ 
1 + u + 2yu' ' ’ 


1 + U 


+ 2u + (2 — d)yu'. 


(33) 


Note that, along the same lines, one may absorb Ad- 

On projecting, we made just one approximation, that of Eq. |^. Before clos¬ 
ing the section, we briefly analyze how, in the N —>■ oo limit, this approximation 
becomes exact. Eq. p8| amounts to. 


V2{k) 


\k^ + V2{0) 

Um(0, ...,0) m > 2. 


(34) 


Hence, it becomes exact when couplings do not depend on k (except for V 2 that 
has a term). If one expands Eq. in powers of momenta, one realizes that 
an hypothesis of such a kind does not hold for N finite, but there are some 
chances at N = oo. 

To proceed, it is necessary to rewrite Eq. introducing 0{N) invariants, 

(35) 

Then, Eq. ^ at = oo becomes closed for the invariants of the form = 

dki+k2 : 

S = hm — lii{2-—)+dS+(2-d) / ^kj^- / ^kk- — — —hconst, (36) 

r^o2TJk S^k Jk d^k Jk dkS^k 














with S is the piece of S that contains only the invariants (called diagonal 

piece in Ref. Q), and we set 77 = 0 which is a general result on the N —f 00 
limit. 

An inspection of how couplings of the diagonal piece at k = 0 are affected by 
those with increasing powers of k suggest that, if all couplings at a certain power 
of k are set to zero, all the rest but the ones at A: = 0 can be simultaneously set 
to zero, so it should not come as a surprise that an ansatz like 




(37) 


with holds. 

Indeed, Eq. ^ leads to a closed equation for V{(p^), 

V = Ad In (1 + + dV + (2 — d)ip^V '+ const, (38) 

which turns out to be the same as the large N limit of Eq. |^. This is the result 
we were after. The LPA in the large N limit is not only solvable as we will see 
in Section |^, but also exact. If our initial potential has couplings that do not 
depend on fc, the exact evolution for those is governed by Eq. 38. 


4 Results 


Fixed Points (FPs) u*{y) of Eq. ^are the analytical solutions of 



2>u' + 2yu" 

1 + u + 2yu' 


+ (IV-1) 


1 + M 


+ 2u + (2 — d)yu', 


(39) 


with Ad some d-dependent constant (Eq. |3^) which equals (27r^)“^ in d = 3. All 
universal quantities are independent of its precise value. Note that this is an 
ordinary differential equation of second order and, thus, two initial conditions 
are required. Nevertheless, at y = 0 Eq. ^ is singular, since the coefficient of 
its highest derivative vanishes, and we must supply the analyticity condition 


( 0 ) = - 


2N 


AdiN + 2) 


uUl + K), < = ^*(0), 


(40) 


which reduces the possible solutions to a one parameter set. 

A careful study shows, however, that for any initial value Uq — the solu¬ 
tion remains negative for any y, resulting in an unbounded potential. We, thus, 
end up with Uq > — 1 as only possible values. Furthermore, both analytical and 
numerical considerations show that solutions of Eq. ^ end up with a singularity 


for all but a finite number of Uq. In Fig. we plot the singular point j/s as a 
function of the initial condition Uq for different values of A^ at d = 3. For d < 4 
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8.0 
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0.0 ^ = 

- 1.0 2.0 5.0 8.0 


u^(0) 


Figure 3: Point ys where the singularity is encountered as a function of the FP 
initial condition u*(0), for fV = 1 (dashed line) and N = A (solid line), in d = 3. 


only two initial values are allowed, one corresponding to the Gaussian Fixed 
Point (GFP) and another one to the Heisenberg Fixed Point (HFP)^ We study 
each in turn. 

The GFP is 

u*{y) = 0. (42) 

Its eigenoperators are solutions of the linearized version of Eq. ^ around rt = 0, 

^9 = + {N — l)^'] + (2 — d)yg' + 2g, (43) 

with g{y) the derivative of the scaling operator and A its eigenvalue. For the 
sake of convenience, we define 


z = 


{d-2)Ny 

2Ad 


(44) 


^In fact, there is a third one, Uq — s- +oo. It corresponds to the Trivial Fixed Point with 
zero correlation length (infinite mass). 
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to obtain 


zg'\z) + (1 + y - z)g'{z) + J^giz) = 0, 


(45) 


which is the generalized Laguerre equation. Polynomial solutions select A to be 

A = 2-n(d-2) (46) 

with n a non-negative integer, and, in such a case, 

g{y) = LW^H^y) (47) 

with the generalized Laguerre polynomial |^. Critical exponents are 

the Gaussian ones, which allow to define massive interacting QFTs in the vicinity 
of this FP. 

The special case of = 1 was previously studied in Ref. 0 in terms of the 
variable x = The solutions are the Hermite polynomials of odd degree, 
H 2 n+iiVz). Indeed B, 


H2n+li\/z)- 


On the other hand, for N 


n!22"+iy^ 
oo the eigenvectors simplify to 
1 / A, " 


»'*'> = ;n I" 


d-2 


o{i), 


(48) 


(49) 


a result which will be derived independently in Section p[ A direct diagrammatic 
computation of the above results is left to Appendix 

The HFP has to be studied numerically. From the asymptotic behaviour of 
Eq. ra, it is, for y —> oo, 


I'iy) 


B- 


2A, 


-.y 




d{d-2y 
2Ad (. 4 


(50) 


(d2 - 4) 


1 - 


{dZr2)N 


y 


d+2 
■ d-2 


with B an a priori unknown numerical constant. Imposing this dependency 
for large enough values of y, together with the consistency condition at ?/ = 0, 
selects the true FP out of the whole bunch of non-analytical solutions. 

Recall, however, that Eq. ^ is singular at the origin. This makes it difficult 
to numerically integrate the equation from large y towards the origin, although 
it is quite easy to integrate it from the origin. We circumvent this difficulty by 
shooting to a fitting point |10 . That is, we integrate forward to some point yo 
from a tiny value of y, satisfying the condition in Eq. M- Which tiny value to 
take is immaterial as long as one checks it is less than the allowed errors. At 
yo we compare it with the backward integration from a large value where the 
asymptotic condition is imposed. In this way, the constants Uq and B from Eqs. 
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3 ^ 


Figure 4: Solutions corresponding to the Heisenberg Fixed Point in d = 3. The 
long-dashed line corresponds to iV = 1, the dashed one to iV = 4, the dotted 
one to IV = 100 and the solid one to N = oo. 


010 are fixed. The precise value of i/q, of course, does not matter, but we 
found it efficient to take values close to ^ 0.05, where the potential is known to 
have a minimum in the exact solution at iV —> oo. Plots for some values of N 
at d = 3 are shown in Fig. where they are compared with the corresponding 
large N solution. 

Linearizing Eq. around the HFP, its eigenvectors are found to fulfil 


^5 


Ad 


3g' + 2yg" 
l + u* + 2yu*' 

9' 


(3 u*' -b 2y u*"){g -b 2yg') 


(1 -b u* -b 2y u*'Y 


+ (iV-l) 


1 -b u* 


-{N-l) 


u*' g 
(1 -b 


+ (2 - d)yg' + 2g, 


(51) 


which can be 
exponentially, 


solved using similar techniques. In general, solutions will 


grow 

(52) 
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N 

V 

UJ 

N 

V 

UJ 

N 

V 

UJ 

0 

0.6066 

0.5432 

7 

0.9224 

0.8876 

50 

0.9895 

0.9861 

1 

0.6895 

0.5952 

8 

0.9323 

0.9028 

60 

0.9913 

0.9884 

2 

0.7678 

0.6732 

9 

0.9401 

0.9145 

70 

0.9925 

0.9901 

3 

0.8259 

0.7458 

10 

0.9462 

0.9238 

80 

0.9935 

0.9941 

4 

0.8648 

0.8007 

20 

0.9736 

0.9639 

90 

0.9942 

0.9924 

5 

0.8910 

0.8396 

30 

0.9825 

0.9764 

100 

0.9948 

0.9931 

6 

0.9092 

0.8673 

40 

0.9869 

0.9825 

1000 

0.9995 

0.9994 


Table 1: The critical exponents v and lo for different values of N and d = 3. All 
digits are significant. The N ^ oo result is :/ = a; = 1. 


N 

A 2 

A 3 

N 

A 2 

N 

A 2 

N 

A 2 

1 

-2.8384 

-5.1842 

5 

-2.8873 

10 

-2.9426 

60 

-2.9917 

2 

-2.8348 

-5.1100 

6 

-2.9036 

20 

-2.9731 

70 

-2.9928 

3 

-2.8482 

-5.0582 

7 

-2.9167 

30 

-2.9827 

80 

-2.9938 

4 

-2.8680 

-5.0286 

8 

-2.9273 

40 

-2.9873 

90 

-2.9945 


9 

-2.9357 

50 

-2.9899 

100 

-2.9951 


Table 2: The third and fourth eigenvalues, A 2 , A 3 , for different Ns at d = 3. All 
digits are significant. The N ^ 00 result is A 2 = —3, A 3 = —5. 


whereas for special values of A they are much smoother. 


9{y) ~ y'^ + 

(A + 2)Ad 


AAj 


d{d-2)B 


1 - 


2A 


{d? - 4)^2 


1 - 


A^(d + 2) 
4(d + 2) + 2A(d-2) 


d 

y 


(53) 


A^(d + 2)2 


d+2 


This asymptotic behaviour, together with an arbitrary normalization of eigen¬ 
vectors, leads to the quantization of A. The results for v and oj are summarized 
in Table |^, where v = I/Aq and uj = — Ai, Aq, Ai being the first (the relevant) 
eigenvalue and the second (the first irrelevant) one, respectively. With nearly 
the same ease next irrelevant operators can also be studied. In Table we 
present our results. 

The A^ = 1 case was first computed in Ref. 0, with results v = 0.687(1) and 
u = 0.595(1), and afterwards in Ref. |]ll|. Our numbers are compatible with 
the former and coincide with the latter. 

Eigenoperators are obtained as an additional bonus. As a matter of example, 
we plot in Fig. ^ the first eigenvector, normalized to ( 7 ( 0 ) = 1/2, corresponding 
to different values of A^ at d = 3. 

The A^ —!■ 0 limit deserves a special remark, as the RG equation 33 is singular 
at A^ = 0. Nevertheless, one may compute its eigenvalues in an indirect way. 
Taking values of N between 0.1 and 0.01, one finds that they perfectly fit a 
straight line, with no sign at all of subleading corrections. Critical exponents 
reported in Table ^ are just their extrapolation to A^ = 0. In Fig. ^ we plot the 
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Figure 5: The first eigenoperator, normalized to (jf(0) = 1/2, for = 1 (long- 
dashed line), TV = 4 (dashed line), N = 100 (dotted line), N = oo (solid line), 
at d = 3. 

first eigenvalue for 0.01 < N <0.1 together with its linear extrapolation. 

This apparent singularity at iV = 0 is, however, due to the rescaling in Eq. 
Had we defined the partition function without explicit factors of N, the 
Ad/N coefficient in RG Eq. ^ would be simply Ad, allowing a direct analysis 
of the N < 0 values. We just decided to stick to our normalization. 

5 Comparison with other equations 

In this section we comment on other ERG equations and compare with known 
results. 

Apart from the Wegner-Houghton equation studied above, one of the oldest 
ERG equations appeared in Ref. |^. It was further modified by Polchinski 
to put it in the form we now use. The regulator is chosen in such a way that 
every propagator is multiplied by a function K{p^) which is analytic everywhere 
in the finite complex plane, vanishes faster than any polynomial for large real 
values of its argument and it is standardly normalized to K{Q) = 1 ||^. The 
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Figure 6: The first eigenvalue for small values of at d = 3. The dots are the 
calculated numbers, the straight line is fitted to extrapolate to fV = 0. 


evolution of the action takes the form 


S 


5S . SS 1 
(50 ■ ■ d0 “ iV 


+ dS + 


2 — d — ri{t) 
2 


K'^ 


S^S 

S(l)5(j) 



/ SS 


(54) 


where a dot means both summation over discrete variables and integration over 
continuous ones. As the non-linearities are limited to quadratic terms, the 
projected equation presents a nice form quite amenable for numerical analysis. 


+ (^1 +^ + {2-d)y - 2yu^ u' + {2 - u)i 


(55) 


where a convenient rescaling of the potential and the field variable y has been 
made. It has been previously studied for = 1, using the variable x = in 


Ref. [IR 


The GFP is also u = 0, with its eigenoperators satisfying Eq. M after the 


15 






rescaling z = (d — 2)y. The solution u = 2 corresponds to the Trivial Fixed 
Point, ^ also known as High Temperature Fixed Point. 

The analysis of the HFP is analogous to the Wegner-Houghton one. There 
are singular solutions which behave 


2 1 
N Vs-y' 


(56) 


whereas the true FP must fulfil the conditions 


u*{y) 


- 2-By~^ 

y—too 


(57) 


with B a non-vanishing constant. Again we can shoot to a fitting point to obtain 
the FP. Eigenvalues A are obtained from the linearized form of Eq. imposing 
the eigenvectors to be y-i>'+ 2 )/(d+ 2 ) y Some of them are shown 

in Tables Si We leave to Appendix ^ the study of its large N limit. 

Another quite used equation is based on the evolution of an infrared-regulated 
one-particle-irreducible effective action [||, 


=4“' 


K- 


dK 

'M 


(^1 + K- 



(58) 


where K is some regulating function and r[(p] is the effective action minus the 
kinetic term ^(p-K~^ -ip. It presents some attractive features: in the limit where 
the IR cutoff vanishes (which corresponds to the continuum limit) it is directly 
observable, as it is the equation of state of the system. Moreover, it presents 
a natural expansion which can be directly mapped to the usual loop expansion 
of perturbation theory ]l^ . Various regulators have been used in connection 
with this equation. If a sharp momentum cutoff is chosen, its LPA coincides 
with that of Wegner-Houghton’s. On the other hand, it has been extensively 
worked out using the special polynomial cutoff introduced in Ref. [^, which 
has the advantage of preserving some symmetries of the exact equation through 
the derivative expansion. This equation is much harder to numerically analyze 
than the above two ones. We do not discuss it further, but only quote the results 
(Tables ^ and refer the reader to the original literature ©H- 

To assess the systematic errors involved in such calculations would require 
to compute corrections to the LPA. Instead, we compare our results with some 
derived from other techniques. 

There are accurate determinations for v from a variety of methods such as 
£ expansion ||l^, Monte Carlo Renormalization Group ||l^, |^, Monte Carlo 

“^To substantiate it, consider a potential consisting only of a mass term, u = Note 

that at first order in r (or, alternatively, for t —> —oo) it coincides with the relevant mass 
operator from the GFP, u = 0, whereas ^ for t —> oo. Furthermore, it is easily 

proved that all linear deviations from u = 2 are irrelevant. 
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(3) 


N (1) (2) 


(4) 


0 

0.6066 


0.5880(15) |17| 

1 

0.6895 

0.6496 

0.6604 

15| 

0.624(1)(2)]^ 

2 

0.7678 

0.7082 

0.73 1 

16 


0.6721(13) 1 

20 


3 

0.8259 

0.7611 

0.78 I 

16 


0.7128(14) 

20 


4 

0.8648 

0.8043 

0.824 |161 

0.7525(15) 

20 



Table 3: Critical exponent v, &t d = 3, for different values of N, using different 
RG equations (all digits are significant): (1) Wegner-Houghton, (2) Polchinski, 
(3) effective action; together with accurate determinations (4). 


N (1) (2) (3) _^ 


0 

0.5432 


0.80(4) 1 

17 


1 

0.5952 

0.6557 

0.6285J 

15| 

0.85(5) I 

18 


2 

0.6732 

0.6712 

0.66 1 

16 

1 

0.780(25) 

17 


3 

0.7458 

0.6998 

0.71 

16 

1 

0.800(25) I 

17 


4 

0.8007 

0.7338 

0.75 

16 

1 

? 


Table 4: Critical exponent w, at d = 3, for different values of N^ using different 
RG equations (all digits are significant): (1) Wegner-Houghton, (2) Polchinski, 
(3) effective action; together with existing estimates (4). 


simulations and Strong-Coupling series |^. Unfortunately, estimates 

for oj are not so precise. A comparison of our results from Wegner-Houghton 
and Polchinski equations with those from the IR-cutoff effective action and other 
existing estimates are shown in Tables 

Critical indices from Wegner-Houghton equation are 10—15% off for N < 10, 
except at = 0, which is quite accurate (^ 3%). In the interval 10 < iV < 100, 
they are off by a few per cent, and for N > 100, the errors are less than one 
per mile. Other ERG perform slightly better, specially Polchinski’s equation in 
which critical indices are a few per cent off at most. This should not come as 
a surprise, since in spite of its conceptual clarity, the sharp momentum cut-off 
leads to bad locality properties of the renormalized actions, being then very 
much sensitive to truncations. In order for the projected ERG equation to 
deliver more precise outcomes, one should find RG blockings with the property 
that renormalized actions are as much local as possible, as it is achieved in the 
lattice regularization with blockings in real space . 

6 N ^ oo 

Wegner and Houghton |Q already studied the limit TV ^ oo of their equation, 
encountering both the GFP and the HEP. They also briefly mentioned the possi¬ 
bility of further FPs. Our aim is to discuss in certain detail those not previously 
studied and present a general framework to deal with this kind of large N equa¬ 
tions. An application of these methods for the Polchinski equation is given in 


17 






































Appendix g. 

The projected RG equation becomes first order in the N ^ oo limit, 


u = Ad 


1 + M 


+ (2 — d)yu + 2u, 


which is easily solved for 2 < d < 4 once we invert u and consider y{u, t), 

Ad 


y = -2uy' - (2 - d)y - 


1 + u' 


( 59 ) 


(60) 


The prime now stands for derivative with respect to u. Its general solution is 
y{u, t) = + f{u), (61) 

with an arbitrary function h{z) and 


f{u) = Ad 


,3-d 


d-2 


+ u dz 


l + uz^ 


(62) 


which is analytic for rt > — 1. Although derived differently, this is the form of 
the ERG equation used in Refs. ^ for d = 3. 

Returning to Eq. M, FPs verify 


2u 


A, > 


(d-2 )y-^ 


(63) 


which can be singular when the denominator of the rhs vanishes. But, if this is 
the case, then analyticity of the FP would imply m = 0, and, therefore, these 
would-be singular points reduce to a unique one. 


y = 


d-2 


, u{y) = 0. 


Furthermore, we can Taylor expand the FP around ji, 

OO 

u{y) = a(y - d) + X! ~ 


(64) 


(65) 


k=2 


an expansion that will be used throughout. 
Eigenoperators are 


with 


ix{y) = 


Adu' 


+ X-2 


dz tx{z), 

(66) 

, (d 2)y 

1 + u 

(67) 


Eq. ^ fixes the only singularity of t\{y) to be a pole at y, with residue 


1 


A-2 


1 + ay \2 — d 


ay 


( 68 ) 
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For eigenoperators to be single-valued, the residue must be a nonnegative integer 
n, which fixes 

An = 2 -I- (2 — d) [ay -I- n(l -I- ay)]. (69) 

Note how critical properties of the system are completely determined by the 
behaviour of the FP around this special point y. 

After those general considerations let us study the different FPs in turn. 
The easiest one is the GFP, which corresponds to 

'^{y) = 0=>a = 0=>A„ = 2 — n(d — 2) (70) 

and, from Eq. ^ 

9{y) = {y-vT, (71) 

a result previously obtained in Eq. 

The HEP corresponds to h{z) = 0 in Eq. 

y(u) = /(„) = A, ^ 

In this case, a = Hence, 

A„ = d — 2 — 2n, (73) 


which are the critical exponents of the spherical model. Eigenvectors are ob¬ 
tained plugging Eq. in Eq. Note that in the limit d ^ 4, the HFP merges 
with the GFP, as can be seen both from Eqs. YA and 73, 


However, contrary to the finite N case, there are more FP. They correspond 
to u) becoming t-independent. 


a — ^ 

y{u) - y = Cu 2 -I- AdU 


^3-d 


dz- 


1 -I- uz'^ ’ 

with C a, so far, arbitrary constant. The behaviour around y. 


u(y) = 


y-y 

C 


+ o((y-y)A), 


(74) 


(75) 


defines for 2/(d — 2) = n, n = 2,3,... a line of FP labelled by the parameter 
C. Since a = 0, critical indices coincide with those of the GFP, for all allowed 
values of C. 

Before getting into a detailed analysis, let us state that, for any allowed C, 
these FPs are reached by Eq. ^ at t > oo, from an initial potential 


V{y) = 


C 

2(72 + 1 ) 


^ y - nAd/2 '\j"^^ 


(76) 


which shows that those FPs encode the large distance properties of local actions. 
Note that C = -too corresponds to the GFP. 
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To proceed further, we should distinguish between n odd and n even. The 
FP in the n odd case is 


y = 


Cui + Ad 



dt 


1 + uf^ ) 


(77) 


The admissible values of C are 0 < C < oo. Although the HFP corresponds 
to C = 0, it does not belong to this line of FPs. This is not so surprising in 
view of Eq. critical properties are not continuous in C at C = 0. Note, 
however, that these singularities have nothing to do with RG transformations: 
C labels different initial conditions for the RG transformations, but once one C 
is chosen, it is not changed along the flow. 

The FP in the n even case comes from matching the two branches of y{u) 
at y. 


f -Cui + f{u) 0 <y<y 
\ TCu" + f{u) y>y 


(78) 


Now the line of FPs extend from C = +oo up to the end point 


C = 


AdTT 
2 sin 21 ’ 


(79) 


which appears as the existence condition of the FP at y = 0. 

The case C = deserves further considerations because the FP is not 

analytical at y = 0. It behaves like u{y) = y +^(1)> eigenvectors 

like, 

(n + l)(n + l-k) 

gk{y)^y • (80) 

The case n = 2 (d = 3) were previously studied in Ref. and these features 
were identified as the RG structures which make the existence of the BMB 
phenomenon possible. That is, the possibility that at C = the FP 

acts as an ultraviolet stable one, allowing to define a new continuum limit other 
than the Gaussian and Heisenberg FPs. In fact one that has mass generation 
and spontaneous symmetry breaking of scale invariance accompanied with the 
appearance of a dilaton. What we now see, is that all RG properties that make 
this phenomenon to appear at d = 3, are present for all the rest of n even values, 
hence a similar scenario is likely to take place. To elucidate these issues would 
require a further analysis of the effective potential at the vicinity of this peculiar 
FP. This is, however, beyond the scope of this paper. 

The preceding considerations have shown the existence of a line of FPs start¬ 
ing at the GFP and including it. This scenario is only feasible provided that the 
GFP contains marginal non-redundant operators. This restricts the dimensions 
of space (space-time) to be 


d = 2 


n = 2,3,... 


(81) 


which gives a nice physical interpretation of the dimensions previously derived 
from analyticity considerations. Thus, in the large N limit, the marginal opera¬ 
tor becomes completely marginal and a line of inequivalent FP appears, though 
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they have the critical exponents of the GFP. This might be a general feature of 
this limit. 

Before closing this section we would like to remark that, although our anal¬ 
ysis of Eq. 60 is restricted to dimensions 2 < d < 4, a similar study in other 
dimensions is possible along the same lines. For instance, in d = 4, Eq. ^ is 
substituted by 


2t 


y{u, t) = 


u 1 

T ~T 


Stt^ IOtt^ 


dz 


1 -I- uz"^ ’ 


(82) 


which does not have any FP other than the GFP, u = 0. On the contrary, at 
d = 2, 

y{u,t) = h{e~'^*u) - ^ + (83) 

ZTT AtT 

where a FP appears for h{z) = —^\n{z). A thorough study for all d would 
take us too far afield. 


7 Conclusions 

This paper presented a detailed study of Wegner-Houghton equation. We have 
shown that in spite of its intractable appearance, the use of the LPA makes it 
amenable for extracting plenty of valuable information. 

The derivation of the equation has been worked out, with special care about 
the terms that contribute to the blocking. The projection has been derived, and 
explained in which sense it is exact in the N oo limit. 

What we feel remarkable about the ERG methods within the LPA is, suc¬ 
cinctly; (i) The qualitative features of the RG at 2 < d < 4 are the same as in 
the expected full computation; (ii) Eigenvalues and eigenvectors are obtained 
with ease, both relevant and irrelevant ones; (iii) Numerical results are, at this 
order in the derivative expansion, not as accurate as other methods, at least for 
low enough N, but those are computationally much harder; (iv) The N ^ 0 
limit is reliable within this framework. 

The LPA of Polchinski equation has been also studied, and found that the 
accuracy is improved. This suggests that a calculation of next order in the mo¬ 
mentum expansion may deliver quite rigorous results. Nevertheless, ambiguities 
appear in connection with the determination of the non-vanishing anomalous 
dimension, which are not present at lowest order fl^ . A further study of these 
issues, hence, seems justified ^T\ . 

Section ^ has been devoted to the N ^ oo limit. Some existing results, scat¬ 
tered throughout the literature, have been reviewed within a general framework, 
which has also been fruitful to study new FPs at critical dimensions d = 2-|-2/n, 
as well as to deal with other RG equations. As an example, we have worked out 
the large N limit of Polchinski equation. 

The scope of the paper did not permit to consider other subjects of interest. 
Among those, the study of the projected RG equation outside the linear 
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approximation, useful to elucidate the topology of the RG |^, ^ ; the study of 
the effective potential of the new FPs at > oo, with special emphasis on the 
singularities that we associate to the BMB phenomenon; and a more systematic 
study of Polchinski equation. We hope that those issues will be tackled in a 
near future. 

Before closing, we would like to come back to our original motivation on first- 
order phase transitions, singling out some of our results which have implications 
for them. We have explicitly shown in Section ^ that in the large N limit, the 
LPA is exact (see Eq. ^). This equation, in the form of Eq. was shown in 
Ref. 0 to lead to singularities of the RG whenever two competing minima ex¬ 
ist, which is the case when first-order phase transitions take place. Those exact 
results point that the RG trajectories are singular and multi-valued, which is at 
odds with Fundamental Theorems One and Two of Refs. [||. Nonetheless, one 
should keep in mind that those theorems do not rigorously apply for unbounded 
spins. On the other hand, the derivation of Wegner-Houghton equation of Sec¬ 
tion || assume positivity of the eigenvalues of at the step of Eq. |l^, which, 

as suggested in Ref. Q, may not be the case if metastable states do exist. A 
further study on all those subjects is beyond the scope of this paper, but we 
think that they pose interesting questions on the conventional wisdom of the 
RG. 

Acknowledgments 

It is a pleasure to acknowledge interest and discussions with D. Espriu, J.I. 
Latorre, S. Leseduarte, T.R. Morris and F. Sala. The contribution of J. Herrero 
and A. Pineda in our discussions is also acknowledged. J. C. acknowledges H. 
Osborn and DAMTP, where this work was initiated, for kind hospitality. A. T. 
has benefited from a fellowship from Generalitat de Catalunya. This work has 
been supported by grants AEN95-0590 (CICYT) and GRQ93-1047 (CIRIT). 


After the completion of this paper, we became aware of the work of Ref. , 
where a computation of the exponent z/, for different values of N, is performed, 
using the LPA with a further truncation in the number of fields. Although this 
additional expansion is known to be problematic |p^ , their results are in good 
agreement with ours. 

A recent reference about the large N limit and the BMB phenomenon is . 
We thank M. Moshe for pointing it to us. 
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Figure 7: Diagrams to be considered, where solid lines mean integrated modes 
and dashed lines are not. A, B, C, are the contributions of the quartic vertex 
to the identity, the mass operator and the quartic coupling itself. D is a generic 
diagram with L loops. 


A The Gaussian RG Equation 

In this section we work out the GFP from a direct computation of the partition 
function, reobtaining in quite a different way expressions already quoted in 
Section We use diagrammatic techniques very similar to the ones in deriving 
the approximate recursion formula |^. We explicitly treat the = 1 case, but 
the N arbitrary situation is completely analogous. 

We consider the splitting in Eq. 

but now accounts for modes with i < 9 < 1- Then, the dilatation in Eq. 

0is 

k = 2q, = (85) 

We consider our action as the GFP and the subspace spanned by (ultra-)local 
operators. 


S[(t>]=1+ \ j (tiqq^(t>-q+ X! 9^''^ [ (86) 

dq 1=2,4,■■■ 

where we explicitly include the identity. 

The linear approximation amounts to consider g^') couplings to be infinitesi¬ 
mal. In such a case we can compute perturbatively the flow. In Fig. |^, the three 
diagrams (A,B,C) generated from a quartic coupling are shown. A generic dia¬ 
gram of a 2n operator with L loops, like the one in Fig. ^ D, gives a contribution 
to the 2n — 2L coupling 


(2n)!c~^ o2ri-2L+d(l-ra+L) 

{2n-2L)\L\2^ 


(87) 
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where c = The linear RG matrix is upper triangular, with non¬ 

vanishing matrix elements, 




(2fc)!c^ ^ 

( 20 !(fc- 0 ! 


( 88 ) 


The main result is that the eigenoperators of the linearized RG matrix are 


On(x) = 

1=0 


{2n)l 


^(20!(n-Z)! 


2(22-d _ 1) 


n—l 


•p (x), (89) 


with eigenvalues A„ = 2'*'" = 2^^ n)d+ 2 n_ proof is by direct substitution, 


„(fe) _ o3l-n+d(l-l) {n-l) i‘^ri)\ ^ (2^ _ l)fc 


_ ty2n+d{l — n) ' 

(20!(n-0! V2(22-‘i-l) 


= A„5«, (90) 


(n- 0 ! 


where the relation 

n 

^^jk-iy.{n-ky. 

has been used. Eq. ^ may be now rewritten as 

f Ad 


i,_n \ / 




/A 


(91) 


(92) 


since c = / -^ = 7|^(1 ~ 2 *^^ d)y g^j.g |.]^g Hermite polynomials, 

which, after using the relation ^ 

H2n{x) ='^nH2n-i{x), (93) 

reduce to the result already quoted in Eq. |4|. 


B N oo with Polchinski equation 

In this appendix, we sketch the results for the large N limit of Polchinski equa¬ 
tion. Our aim is twofold: on one hand to show that our techniques introduced 
for Wegner-Houghton equation are quite general and may be easily applied to 
other RG equations as well; also to prove that the annoying restriction it > — 1 
of the former equation is not an essential feature. This last fact would make 
interesting to repeat the large N analysis of Ref. 1^ for Polchinski equation. 

At N ^ oo, Eq. ^ becomes 

It = [1 — (d — 2)y — 2yu\u' + (2 — u)u. (94) 
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Its general solution is 



( 95 ) 


with h{z) an arbitrary function and 



I 



d-2 

2 


dz 


(96) 


Note that f{u) is analytic for0 u < 2, unlike the Wegner-Houghton case, in 


which analyticity has a lower bound at u = — 1. 

For u = 0 the point y = satisfies u(y) = 0. Expanding the FP solution 
around it, u(j/) = a(y — y) + O ((y — y)^), and imposing the eigenoperators to 
be analytic at y, y(y) ^ (y — y)", the condition for eigenvalues equivalent to Eq. 



IS 


An = 2(1 — ya) — n{d —2 + 2ya). 

The GFP corresponds to u = a = 0, with eigenoperators 


(97) 



(98) 


The HFP corresponds to h{z) = 0, which gives a = and A = 

d — 2 — 2n. 

There is a line of FPs for d = 2 + 2In, 


d 



(99) 


which starts at (7 = +oo (GFP). It ends, for n odd, at C = 0 (but without 
reaching the HFP) and, if n is even, at C = Cq = Precisely C = Cq, 

for n even, is the only value where the FP is analytical at it = 2 (it is Taylor- 
expandable around y = with = 2). This allows u = O(^) in the 

vicinity of y = 0, like in the Wegner-Houghton case. These non-analyticities 
associated to the BMB phenomenon seem to present, therefore, some kind of 
universality. 


®Recall that u = 2 corresponds here to the Trivial Fixed Point. 
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